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Abstract 

We show estimates for shadowing of gluons at small values of x, appropriate 
to RHIC and LHC experiments. Using a new evolution equation which takes into 
account the effects of gluon recombination to all orders in gluon density, we show 
that there is a significant depletion in the gluon density of large nuclei. 
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1 Introduction 



Understanding the initial conditions of a relativistic heavy ion colhsion from first principles 
is perhaps the single most challenging problem facing the heavy ion community. Proposed 
signatures of a possible Quark-Gluon Plasma (QGP) formed in heavy ion collisions will 
crucially depend on the initial conditions. At the very early stages of the collision, one 
would need to take the full quantum mechanical nature of the nuclei into account which 
is a prohibitingly difficult task since it would require full knowledge of the nuclear wave 
functions. The McLerran-Venugopalan model offers a new and promising tool to inves- 
tigate these early times by representing the initial nuclei by classical fields. At later times, 
when the highly off-shell modes of the field are freed by hard scattering and go on-shell, 
one can identify these modes with partons and consider the initial distributions of these 
partons in the nuclei. 

When calculating a typical nuclear cross section, one needs to know the distribution 
of a given parton kind in the nucleus jgjA- Using the QCD factorization theorems, one can 
then write the nuclear cross section as a convolution of the nuclear parton densities with 
the hard parton-parton cross section, i.e. 



In order to obtain the nuclear parton distributions, one can take the corresponding 
parton distributions in a nucleon and scale them by the atomic weight A. This sounds 
plausible specially at high values of since one does not expect nuclear effects to be 
important at large values of . However, this expectation was proved to be too naive and 
it was experimentally found the the distribution of partons in free nucleons are strikingly 
different from those in bound nucleons. This difference is more pronounced in large nuclei 
at small values of x^j where a significant depletion in the number of partons is observed so 
that a simple A scaling does not hold. 

Alternatively, one could measure the nuclear parton distribution in a particular exper- 
iment. Since these distributions are universal, one could then use them to predict nuclear 
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cross sections in other experiments. The most recent experiments measuring nuclear par- 
ton distributions have been performed by the NMC collaboration at CERN SPS and 
by the E665 collaboration at Fermilab 0. However, both of these are fixed target ex- 
periments and are limited in the kinematic range in x and they can cover. Also, the 
amount of data in the kinematic region where perturbative QCD would apply is limited. 
A lepton-nucleus collider would go a long way towards expanding our knowledge of the 
nuclear parton distributions and is urgently needed. 

Once the nuclear parton distributions are known at a given value xq and Qq, one 
can use the perturbative QCD evolution equations to predict the distributions at different 
X and Q. However, the standard evolution equations are expected to break down at 
very small values of x due to parton recombination effects. A new evolution equation 
(JKLW) which takes these effects into account was derived in 0] and is the non-linear all 
twist generalization of the standard perturbative QCD evolution equations such as DLA 
DGLAP and GLR/MQ (see also ||^ for a similar equation). These non-linear effects were 
investigated numerically in and were found to be important in the kinematic region to 
be explored by the upcoming experiments at RHIC and LHC. One can also use this new 
evolution equation to predict the x, Q, bt and A dependence of the gluon shadowing ratio 
defined as 



where xGa and xGn are the nuclear and nucleon gluon distribution functions respectively. 

2 Shadowing of gluons 

To calculate the shadowing ratio S for gluons, we start with the following evolution equation 
which was derived from the effective action for QCD at small x in [Q. 
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and Ei(x) is the exponential integral function defined as |^ 

poo g-(l+i)2: 

Ei(s) = / dt , X > (5) 

^ ' Jq l + t ^ ' 

This equation was shown to reduce to DLA DGLAP and GLR/MQ at the low gluon 
density limit. In we showed in detail how to solve this equation numerically. Here, we 
briefly review our main approximations and assumptions and refer the interested reader 
to for more details. In order to solve equation (|^), we need to know the initial gluon 
distribution at some reference point xq and Qq as well as its derivative (this is due to 
making the semi-classical approximation). We then use a fourth order Runge-Kutta code 
to calculate the gluon distribution at any other point x and Q. In ||^, we took xq = 0.05 and 
Qo = 0.7, but the effective Qo for most points calculated was about IGeV. The reason for 
our choices were two fold; first that experimentally it is known that the shadowing ratio is 
about 1 in the range x = 0.05 — 0.07. Also, in order to maximize the effects of perturbative 
shadowing, we needed to start from as low value of Qo possible while keeping it high 
enough so that perturbative QCD is still valid. With these approximations, we showed in 
0] that for large nuclei, the non-linearities of the evolution equation are very important. 
Our results for the gluon shadowing ratio as defined in @ are shown in Figure 1. It is 
clear that gluon shadowing at RHIC and LHC will be important. 

Shadowing of gluons in nuclei will have significant effects on the measured observables 
in the upcoming experiments at RHIC and LHC. For example, initial minijet and total 
transverse energy production will be greatly reduced. Also, heavy quark production will be 
significantly effected since its cross section is proportional to the squared of gluon density. 
Basically, any production cross section which involves the distribution of gluons in nuclei 
will be modified. Therefore, it is extremely important to understand shadowing more 
thoroughly. For example, the observed shadowing ratio (of F2) does not seem to have 
a significant Q dependence while parton recombination models tend to predict a strong 
Q dependence of the shadowing ratio. This could in principle be due to assuming no 
shadowing at the initial point, i.e. the point where the perturbative evolution starts from. 
To investigate this, one should include initial non-perturbative shadowing at the reference 
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Figure 1: Gluon shadowing as defined in (0) 

point and then evolve the distributions with both the leading twist DGLAP and all twist 
JKLW evolution equations. The difference between the two would be a clear indication of 
importance of higher twist effects in understanding gluon shadowing. 
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